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I. INTRODUCTION 


The method of moments (MM) has proven to be an effective technique 
for the analysis of microstrip antennas [1-5]. Most of the previous 
work has concentrated on analyzing microstrip antennas near first 
resonance. In the method of moments analysis of the microstrip antenna 
the two most important quantities to be evaluated are: 

1. the self or mutual impedance between the microstrip expansion 
and test modes which comprise the elements in the impedance 
matrix 

2. the mutual impedance between the microstrip modes and a 
vertical current filament in the dielectric, representing a 
coaxial feed, and which comprise the elements in the right-hand 
side vector. 

Chapter 2 finds the exact near zone fields of the microstrip modes. 
These fields are in terms of the Sommerfeld integrals. These fields are 
used to determine the mutual impedances, which are also in terms of the 
Sommerfeld integrals. Using the method of Pozar [6], efficient methods 
for the evaluation of the integrals are presented. Finally, entire 
domain expansion and test modes are described which would permit the 
analysis of the microstrip antenna over a broad frequency range. 
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II. PATCH DIPOLES ON A GROUNDED PLANAR DIELECTRIC SLAB 

A. INTRODUCTION 

To analyze the mutual coupling between patch dipole modes on a 
grounded dielectric slab, the electric field due to a patch dipole mode 
is needed. A general solution to the field equations pertaining to a 
microstrip patch on a grounded dielectric slab is presented in Section 

B. The surface current on the patch is introduced only through the 
boundary conditions; this simplifies the calculation, which is performed 
in the Fourier transform domain. The solution is essentially of the 
Green's function type even though the Green's function is not 
constructed explicitly. It is exact in the sense that both the 
dielectric slab and the ground plane are taken into account rigorously. 
Thus, surface waves and coupling to adjacent antenna elements can be 
accurately determined. An efficient evaluation of the Green's function 
is also discussed. A general expression for the mutual impedance 
between two rectangular patch dipoles is presented in Section C. An 
ej wt time dependence is assumed and suppressed throughout this chapter. 

B. THEORY AND GENERAL SOLUTION 

The geometry under consideration is shown in Figure 2.1. The 
grounded dielectric slab is infinite in extent in the x,y directions 
with uniform thickness t. A microstrip patch is printed on the slab at 
the dielectric-air interface with current density J s . since both 
regions 1 and 2 are source free, an arbitrary field that satisfies 
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Maxwell's equations can be constructed from two scalar functions [7]: 

\j> ^ which generates a TM field and ^ which generates a TE field, where 

i = 1 for the region inside the dielectric, and i = 2 for the region 

outside the dielectric. Both scalar functions . and \j>. satisfy the 

mi ei J 

scalar wave equation: 





= 0 


* 


where k- 


p 

e r k^ in region 1 (i = 1) 

p 

k Q in region 2 (i = 2) 


( 2 . 1 ) 



Figure 2.1. Geometry of a microstrip patch printed on a grounded 
dielectric slab. 
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k Q is the free space propagation constant, and e p is the relative 
dielectric constant in region 1. 

a 

Upon multiplying v|> . and ^ by the unit vector z, one can 
construct the fields as [7] 

1 

E i(r) = - vx.(z* ei (rn + v x v x ( z *mi (r )) * ( 2 . 2 ) 


1 

H i ( p ) = Vx ( z W ?) ) + j^T v * 7 x ( z +ei (r ^ ’ (2.3) 

where = e r e 0 and £2 = e 0 . The explicit form of the field components 
can be found in Appendix A. 

Since the structure is infinite in the x-y plane, the four scalar 
wave functions can be represented in terms of their spectral functions 
(or 2-D plane wave expansions) as follows: 

in region 1 


V r) 


1 00 
4ir 2 ^ 


in region 2 


'W k x’V cosk zl( z+t ) 

_WvV sink zl (z+t) _ 


-j(k x+k y) 

e y dk x dk y > 


(2.4) 


V r) 

_*e2 (? L 


4ir 2 


// 


WvV 

_*el^ k x ,k y L 


- j k * r 

e dk dk , 
x y 


(2.5) 
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where 


k zl = / e r kg " k x ' k y . Re k zl > 0 , Im k zl < 0 , 

k z2 = / - k^ - k^ , Re k z2 > 0 , Im k z2 < 0 , 

- * * ( 2 . 6 ) 

k = xk x + yky + zk z2 

and 

AAA 

r = xx + yy + zz 

It should be pointed out that using Equations (2.4,5) in (2.2,3) will 
automatically satisfy: 

a. the boundary condition at the conducting ground plane: 

A — 

z x Ej = 0 at z = -t ; (2.7) 

b. the radiation condition as r = |r| -*• <»; and 

c. the criterion for the integrals in (2.5) to converge as z + ». 

To specify the fields uniquely, the boundary conditions at the 
dielectric-air interface must be satisfied: 

z x (E 2 - Ej) = 0 at z = 0 , (2.8) 

and 

z x (H 2 - Hj) = j $ at z = 0 
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Without loss of generality, the current density is assumed to be Fourier 
transformable and to have both x and y components. Hence, can be 
written as 

J s (x,y) = xJ $x (x,y) + yJ sy (x,y) 



{xJ sx (k x-V + y V k x>V ,e ’ J<k * X+kyy)dk x dk y 


( 2 . 10 ) 


Use of Equations (2.8) and (2.9) will determine completely the 
spectral functions ip ml . ,i|» . (i=l,2) . The algebraic details are carried 
out in Appendix A. Only the field components are listed below: 


region 1 (-“ < x,y < » , -t < z < 0) 


E xi< ? ) - 


4tt 2 


J k x k zl " 

// I ^ e i ue e ’W 
or 


-j(k y x+ky) 

e y sink z ^(z+t)dk x dk y 


* 


( 2 . 11 ) 


-j “ 

E yl^ r ^ = 4 tt 2 H 

J -co 

-j( k x x+ k y) 

e J smk zl (z+t)dk x dk y , (2.12) 


jK y k zl 

k 1)).+ ib . 

x T el cue e y ml 
o r 
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(2.13) 


-j » (kx +k y) -j( k x x+k . / y) 

E zl fr)= ^l! vT^l e cosk 2l (z + t)dk x dk y , 




~ jk x k zl 


" k w+n 


a)y o '"el y v ml 


-j( k x x+k y y) 


cosk zl (z+t)dk x dk 


y ’ 


(2.14) 


v ? > - ^ >1 


" jk y k zi 


+ k v^*r 


u»u "el T x muI 


-j ( k x x+k y y) 


cosk zl (z+t)dk x dk 


(2.15) 


and 


h 21 (F) - 


-J 7, ( k x' fk .y ) 




4.2 ' sink ,(zn)dk dk ; 

-00 K o J 


region 2 (-<*> <x,y<«,0<z<«) 


(2.16) 


E x 2 ^ = 4tt2 // 


k x k z2 

k y*e2 + ”V 
J o 


-j( k x x+k v y +k 7 oz) 
e x y Z4 dk ^ # 


(2.17) 
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i) 








irll \*» 


uj£. 




(2.^ 


fc„ *e2 ^° 


h,^ -CO 




-3 


-v - -r 2 n ® e o 

t ^rl n, - 


V&&V 


-uV'V *' 42 




,0 


H *2 


A <*> 

J- Si 

( V f'l * <v-« 2 -<* 


.jt>/*v + * iZ dv * dV * 

r "4 \ e 

\ - v - sv 


US\1, 




*„V 


* “ 7. _ + v^rf 


-IW*/^**** ’ 

e ^ 


H 




« -^2 Si 

(r^ ' ^ - 50 


°° \ — ^ ~ ^e2 


u>\>, 




i) 




i£*v& S 


-ji.v^v*' 42 a ^i 


h i2 


(r') s ** 


i 

-2 ii ® v 0 
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where 


'el 


U)U 


(k 2 +k 2 )D 
v x y' e 


^ k x^sy ' Vsxl 


(2.23) 


coy sink .t 
° zl 

^e2 = 2 2. tVsv 

(k>k 2 )D Sy 


Vsxl 


(2.24) 


«"V/ 



Je r k z2 


(k 2 +k 2 )D 
v x y ; m 



(2.25) 



k zl sink zl t 

(k 2 +k 2 )D 
' x y' m 


l k Ax + W 


(2.26) 


D e = k 2l cosk zl t + jk z2 sink zl t , (2.27) 

and 

D m * E r k z2 cosk zl t * jk zl sink zl t • (2.28) 

It should be pointed out that there are two dyadic Green's 
functions, 6 1 and G 2 , associated with the grounded dielectric slab. G 1 
corresponds to the case where the field points are inside the substrate 
(region 1), and G 2 corresponds to the case where the field points are 
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outside the substrate (region 2). Only six components of each dyadic 
Green's function are considered since the surface current on a 
microstrip patch has no z-component. These components can be identified 
from Equations (2.11) through (2.13), and (2.17) through (2.19) as 
follows: 




(OjJ 


0 - 00 


// 9 i e (F,k x ,k y ) e Jk-r dk x dk y 


(2.29) 


where 


i 

CVJ 

* 

It 


_ 

A 

A 

k 

II 

X 

X" 

X 

+ yk 


A 

A 

r 1 

= xx 1 

+ yy 


A 

A 

r 

= XX + 

yy + 


a = x,y,z 
3 = x,y 

and g^ are obtained from the electric field components which typically 
have the following form: 


E^IF) = U [9 1 ax (? ' k x' k y )3 ' S x (k x- k y ) + 9^ f - k x 

(2.30) 


It follows from Equation (2.30) that the calculations of either 
electric fields or mutual impedances will invariably involve the 
numerical evaluation of an integral of the form: 
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- F < k x' k y> 

r = II ~T5 dk x dk y • (2.31) 

-oo e m J 

which, however, can be facilitated by changing to polar coordinates k,<j> 
where 

k x = kcos , (2.32) 

ky = ksi n <|> . (2.33) 

Thus, 

2tt F(k,<(>) 

r = / kdk / d<|» . (2.34) 

C k o e m 

The contour C|< for the k integration is shown in Figure 2.2. The branch 
cuts for the branch points k = ±k 0 are defined by the analytic 
properties that 

/ 2 2 

k Q - k < 0 on the entire top Riemann sheet; 

b) Re k z 2 > 0 in the first and third quadrants; and 

c) Re k Z £ < 0 in the second and fourth quadrants. 

However, k = ± k n are not branch points since the integrand is a 

I T 2 

single-valued function of k zl = / e r k Q -k . The branch cuts are also 
shown in Figure 2.2. 
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Im k 



Figure 2.2. Proper contours of integration and branch cuts in the 
complex k plane. 
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The denominator in the integrand of Equation (2.24) defines the 
surface wave modes. These modes are determined by the roots of 

D g = 0 (for the TE waves) , (2.35) 

= 0 (for the TM waves) . (2.36) 

/ 2 2 

Substituting a = k o t/e r -l and a = t / e r k Q -k into Equations (2.35) and 

(2.36) and rearranging terms yields, respecti vely. 



(2.37) 



(2.38) 


It can be demonstrated that if e is real and £ * 1, the roots of 

r r 

Equations (2.37) and (2.38) are real and located inside the segment 

k < k < /e~k . If N ,N represents the number of roots for the TE, TM 
o r o e m r * 

case, respectively, then it can be shown that 
0 , for a < ir/2 

n , for (n— l/2)ir < a < (n+l/2)n, n = 1,2,3,.... , 

and 

= n+1 , for mr < a < (n+l)n, n = 0,1,2 
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It is noted that the dominant TM mode has a zero cutoff frequency. For 
lossy dielectric (Im e r * 0), the roots just move off the real k-axis 
with the form k = k r -jk.., k.. > 0. The exact root locations can be 
determined by using the Newton -Rhapson procedure. 

For efficient numerical evaluation of the integral in Equation 
(2.34), the new contour is adopted by deforming (as shown in 
Figure 2.2). The integration along is performed by computing the 
Cauchy value of the integrals around the surface wave poles. Assuming 
the surface wave poles constitute an ordered set as (p , *=l,2,...,n}, 

X/ 

Equation (2.34) can be written as 


r 


£im 

6+o + 


7 


+ 


7 

P x +6 


+ 


p -6 

+ 7 + 

p n-l +6 


J 

Pn +S 


2tt F( k , tf>) 

kdk / n D d<t> 
o e m 


n 2 tt 

jir l } d<|> Residue 
*=1 o 


F(k,j,) 

DO 
e m 


k=p ( 


(2.39) 


In the case of lossy dielectric, the integrations from p^-6 to 
p^+6, *=1,2,3,... *n, can be evaluated analytically without indenting 
the contour C ot . This is done by using two terms of a Taylor series 
expansion of D e D m about p^, and by taking the value of the numerator 
F ( k , <(>) at k=p^ throughout the interval. 

In actual numerical evaluation of (2-39), it is found [6] that 6 ~ 
0.001k o is adequate, and the infinite integral should be terminated at k 
no less than 150 k 0 for the self impedance calculation. For the mutual 
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impedance calculation, however, this infinite integral tends to converge 
more slowly as the separation between the two dipole modes becomes 
larger, and hence will involve an exhorbitant amount of computer time. 
Special treatment is required to improve the computational efficiency, 
which is dealt with in the next chapter. 

C. MUTUAL COUPLING ANALYSIS 

In order to employ the method to analyze the mutual coupling 
between microstrip patch antennas, one needs to evaluate the elements 
of the impedance matrix and the voltage vector. Assuming the Galerkin 
form of the moment method is chosen, and the basis and weighting 

A A 

functions (modes) are members of the set (D m = xJ mx + yJ my : m=l,2, 

..., N < »} , the mutual impedance between mode J m and mode J n is defined 
as 

z mn = E m (F > ' . (2.40) 

5n 

and the element of the voltage vector is defined as 

v m * ' E m< F l • , (2.41) 

where E m is the electric field excited by mode m, and J. is the 

impressed current source. As presented in Appendix B, the exact 

expressions for z „ and v are given by 
mn m J 
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^mn “ " 4ir 2 


// 


k x k z 2 


( k y*e2 + V^ cos *mn 

<* r\ 


k k , 
y z 2 

- ( k x*e 2 ‘ ^TT V^ sin ^mn 


J nx’ 


( k J*9 + 


k x k z 2 


k y k z 2 


y v e2 r u>e 


*m2) s1n *nm + ( k x*e2 ’ ^T~ V^mn 


-j(k x + k. v ) 

• e x mn y dk dk 

x y 


J'k . < k K> 


m 


„ -jU x x f + k v y f ) 

. 2 j j , $ m1 e y sink„.,t dk„dk. 

4tt^ jj (oe_e_k_ 1 MTU 


-«» o r zl 


v zl u *x u *y 


where 


mtiI 


j e k 

(k x 2 + k y 2 )D m (kx3 “ + kyJmy) 


a)y Q sink zl t 


(UL - Umv) 


* e2= < k x +k >e "' X ’” y y '" X ' 


* C zl S 1 nk zl t 

*™ 2 = ‘ (k 2 *k 2 )D m + kyJ "- y) 



(2.42) 


(2.43) 


(2.44) 


(2.45) 


(2.46) 
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mx 


m 


A 

r 


J m x (x -y> e 


j(k x+k y) 


dxdy 


m 


-x 


m 


(2.47) 


my 


y m x m j(k x+k y) 

/ J my (x,y)e y dxdy 


(2.48) 


“ X m 

m m 


y„ x„ 
.n ,n 


J nx' - r 1 J nx' ««'•»') 


-y„ -x 
J n n 


-j [x' (k cos<|> +k sin*,,,) + y' (-k„sin6 +k cos6 ) 1 
e L x Y mn y Mnrr J ' x T mn y mn dx’dy 1 


(2.49) 

^ny* /" /" J ny'< x '-y'> 

-y„ -x 
J n n 

-j [x ' (k cos* +k .sin* ) + y' (-k sin<f> +k cos*,,,) 1 
e L x r mn y Y mn' J ' x mn y Y mn dx’dy 1 

(2.50) 


(|) is the angle between the x-axis of mode m and the x'-axis of mode n; 
mn 3 

2x m and 2y are the widths of mode m on in the x and y directions, 
m • / m 

respectively, and likewise, 2x n , 2y n are the widths of mode n in the 
x',y' directions (see Figure 2.3). Also, (x^,y^) denotes the 
coordinates of the feed location with respect to the (x,y) axis centered 
on mode m. 
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<*>, the integrands in Equations (2.42) and 


As k = /£* + i/ 2 
K x K y 

(2.43) become highly oscillatory. This will cause numerical difficulty 

in evaluating z and v m accurately. To rectify the situation, it is 

important to note the asymptotic behavior of these integrands as k tends 
to infinity. Defining 


k x k z2 


0)£, 


ID 0 + 

y v e2 

(oe 0 

♦m2 



, k o k y 

k x^my' 

Vmx 

_ k x k zl k z2 

k x^mx~ k y^my 

1 

D e 

. 

— 

— 


— — 

] 


sink 2 jt 

(2.51) 


and 


k y k z2 

T 2 = k x ^e2 “ V 


i oe / 


k o k x 


k x^my“ k y^mx 


k y k zl k z2 
+ D n 


-'m 


k J +k J 
x mx y my 


sink z jt 


x y' 
(2.52) 


one can rewrite Z mn as 


-J 

I 'mn “ 4 tt' 


z - = ^ " t T l(V cos *mn ' V' s1n *J ‘ V 3 ™' s1n *mn + V' 


ny 


-j(k Y x mn + k y ) 

. x mn y^mn' .. .. 

cos4» mn ]} e dk x dk v . 


x y 


(2.53) 


19 



Noting that k zl , k z2 -*• -jk, sink zl t -► e kt and cosk zl t -*■ e kt as k + », 

■j 2 2 

one can approximate, for large k. 


D e - "Jke 


(2.54) 


D m 3 - J ' k 


(e r +l) e 


(2.55) 


and hence. 


n 0 1 

Tl ‘ ' ^ «. r «) 


2 ( £ r + ^) 2 „ 

k o ^ - k x mx “ k x k y^my 


, (2.56) 


C ° k(e r +l) 


2 ( £ r + l) 2 „ „ 

k o ~ 2 - k y ^my -k x k y J mx » ^ 2 * 57 ^ 


Similarly, 


v = 7r // To e 
m 4tr 3 


-j (k x x f * y f ) 


dk x d ky 


where 


j sink z t 

^3 ue o D m ^ k x^mx + k y^my^ 


and, for large k. 


J^o 1 

T3 ”' 2k ° (v « + VV 

k(e r +l) 
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The asymptotic forms of Tp T2 and will be utilized in the next 
chapter to improve the computational efficiency of z mn and v . 

Finally, the proper choice of the dipole modes 3" , m=l,2,...,N, and 
their corresponding transforms are discussed in Appendix C. 
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III. IMPROVING THE COMPUTATIONAL EFFICIENCY 
A. INTRODUCTION 

As mentioned in the last chapter, the integrals for z mn or v m 
converge slowly (especially for large separation between two dipole 
modes, or a dipole mode and the impressed source) when k ■* «. Moreover, 
as k becomes large, the corresponding integrand will get highly 
oscillatory, and will require very small step size in the numerical 
integration scheme. All these translate into large amounts of 
computer-time, and hence high costs. It is, therefore, of great 
interest to improve the computational efficiency of z mn or v m . As 
suggested by Pozar [6], it is useful to construct an integral that have 
two different representations such that the first can be evaluated 
easily, while the second will have the identical asymptotic behavior as 
that of z mn or v m . It is then reasonable to expect the integral for the 
difference of z mn or v m and the constructed integral to converge much 
faster because of their identical asymptotic forms. The total value of 
z mn or v m can be recovered by adding back the constructed integral which 
is evaluated via its first representation. With such manipulation, the 
overall efficiency for calculating z mn or v m can be improved. 

This chapter is devoted to the construction of such integrals. In 
particular. Section B will concentrate on the integral representations 
which exhibit asymptotic behavior identical to those of z mn and v m . The 
alternative representations (which can be evaluated easily) will be 
discussed in Section C. 
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B. MUTUAL COUPLING IN A GROUNDED HOMOGENEOUS MEDIUM 


In this section, an analysis similar to that of Chapter 2 will be 
repeated for a homogeneous dielectric half space D bounded by a 
conducting plane, where u Q and e Q e e are the constitutive parameters. 

Consider a conducting patch S of current density J $ in the x-y 
plane situated at a distance t from the ground plane in D (Figure 3.1). 
The vector potential A due to J $ is given by 

3 s (r') ds' , z>-t (3.1) 

where k e = w/y 0 e 0 e e , 

AAA 

r = xx = yy + zz , 

A A 

r ' = xx' + yy* , 

AAA 

r $ = r-r' = x(x-x') + y(y-y') + zz , r $ = |r g | , 

A A A A 

r i = r-r' + z2t = x(x-x') + y(y-y') + z(z+2t) 


A(r)= 4^ / 
S 


-jk e r s -jk e r. 
e e 

^s ~ ~ 



Homogeneous dielectric half space D 



( IMAGE SOURCE) 


Figure 3.1. Patch current source 3 in homogeneous dielectric half 
space D. 5 
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Employing the identity 


-Jk R 


1 // 1 e 


-JCUx-x') + k (y-y ' ) + k |z-z'|] 


~R = j2tt F7 


dk x dk y > 

(3.2) 


where R = |R| = |x(x-x‘) + y(y- y‘) + z(z-z')| , and 


k z = 


x K y 


, with Imk z < 0 


(3.3) 


one can write 


1 - 

S < r » = 4 ^ W 
S-" 


-J4>i 

e - e 


dk x dk y ds' , z > -t. 


(3.4) 


where 


+ s = < f ) s (f' >k x ,k y ) = k x (x-x') + k y (y-y') + k z |z| , 


(3.5) 


♦i * ♦ i (p > k x ,k ) = k x (x-x') + k (y-y ' ) + k z (z+2t) , 


(3.6) 


It is well-known that the electric field l at reD is given by 


-jn 0 

M -TT 


" 2 1 
K A(F) + 17 VV • R(r) 


(3.7) 


where n Q = / Uo /e 0 and k Q = <i>/u 0 e c 
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Since 


3 9 3 

7 • Mr) = ^ A x (r) + gy A (r) + 37 A z (r) 


1 “ 1 

1 

~ -j+. -H*~ 

e - e 

= 47" 02k z 

1 

C-j. 

X 

S -00 

1 1 

— — 




-ok 


" -J> s 
e - e 


V ?,) 


> dk x dk y ds* 


(3.8) 


one will find 


vv-A(r) = 


* 3 ^ 9 a 3 

x H + y ay + 2 Ji 


v*A(r) 


(3.9) 


1 

4 ? 


/// j2k z 

s— 


2 -j#*- ~ J <f> .• 

-k x (e S -e ) J sx (r) - k x k y 


(e s -e 1 )J sy (r')+ 


+ y 


_ 2 ■J'f’c ~ J4> j 

Vy<* '* > J sx( r ') * k y < e - e >V r,) 


.+ Z 


" k x k z ( sgn ( z) e s - e 1 ^ J sx (r') - k y k z fsgn(z) e 


-j* 


/ V ? ‘> 


dk dk ds' 
* y 
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Next it follows from Equations (3.4), (3.9) and (3.7) that 


E(r) 


-Jio • 1 


fff 


k o 4tt" J s J _l J2 k z £ e 


2 2 

E e k o - k x 


J sx (F '> 


k x k yV ? '> 



■j*s ”3 4>-j \ a 

e - e ) + y 


( 


2 2 
E e k o - k y 


J sy <?,) - k x k y°sx (F '> 


-j* s -j*,' 


- e 


- j 4> - j <(> . 

" z ( k x k z J sx (F,) + k y k z J sy (F,) ^ ( s 9 n(z) e S - e 1 


dk x dk y ds’ 


(3.10) 


Now one can readily formulate the mutual impedance z“ n between 


dipole mode m and mode n as 


z*L = - / E (r+r‘) • J-(r') ds‘ , (3.11) 

mn _ m n 

n 

where E m denotes the electric field in D excited by mode m with current 
m J 

A A 

density 3 (r) e{J (?) = x J + yO : cf1,2 N, which exists 

III Ot CIa Uj 

on the surface patch S m . The definitions of r and r' are obvious from 
Figure 3.2. Observing that 
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E m (r+F' ) = xE mv (F+r') + yE mw (F+F') + zE (F+F 1 ) 
m ' mx my mz 


= x' E m (F+F 1 ) cos* + E m „(F+F' ) sin* + 
mx ' mn my v r mn 


y‘ -E mv (F+F‘ ) sin* + E mw (F+F‘) cos <* + 

J mx v ' r mn my ' r mn 


zE mz < F+F '> 


°n (F '> * xJ mx'< F '> + y' J ny’ (F,) 


and making use of Equation (3.10), one can rewrite Equation (3.11) as 


d _ Jn o * 1 

z mn " 4ir 2 k ^ j2k z e e 
o 


< e e k o - k x) V - k x k y J ~ 


y my 


J nx' cos<1, mn ' J ny' sin<), mn + 


2 2 ~ ~ 

+ u k - k)J -kkj J ,sin* +J , cos* 

v e o y ; my x y mx nx mn ny T mn 


' ' j(k x x mn + k v y mn ) 
e x m y dk dk 

x y 


<v /V /V 


where J , J , J , and J , are defined in Equations (2.47), (2.48), 
in a Hiy nx ny 

(2.49) and (2.50). 
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Letting the impressed source J. be a current filament of 
amplitude I.., normal to the ground plane, and situated at (x^ 
respect to the center of dipole mode m in the z=0 plane, i.e. 

A 

3.j(r) = zI.6(x-x f )S(y-y f ), -t < z < 0 

the mutual impedance between mode m and the impressed source . 
written as 


v m = l E ( F) • J.(r)dv 
m V m ' i ' 


0 

i t ’l E raz (x r z > d z 


-JVi 00 1 

" J 2 k z c e 


Vz J mx * k y k z J »,y 


-j(k x x f + y f r 


0 

/ 

-t 


' jk z -jk 2t - jk r 
e + e 


dzdk dk 
A y 


or 


v d . 


j^o^i ” 1 




-J2k t— 

4ir 2 k ^ 2k z e e 
0 -» 

k J + k J 
x mx y my 


1 - e 1 


-j(k 


x x f + k y y f ) 


dk x dk y 


uniform 
y f ) with 

(3.15) 
i can be 


(3.16) 
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In order to compare the asymptotic forms of the integrands of z, 

and vl with those of z m „ and in Chapter 2, it is convenient to 
m mn m 

express z^ n and in the following forms: 


mn 


mn 


-J 

4tt 


// 


. cos<t> - J„ w iSin<)> 

nx mn ny T mn 


- T 

* r \ 


J^.sin* + 
nx mn 


+ J ny ' cos %n 


"^x x mn + V^ron ^ .... 17 v 

> e J dk x dk y » (3.17) 


d “ jI i 00 

v m = TT~ ft TX e 
m 4 tt 3 

—00 


d -j (k x x f + k y f ) 

0 - x T y T dk.d 


x ky 


(3.18) 


where 


T, » 


-"o 


J Vz £ e 


< e a k o - k x> v - k *yv 


(3.19) 


T d _ n o 
2 " j2k o k z E e 


< £ e k o - k J> - k M. 


(3.20) 


and 


d ~ n o 


T = 

1 o 


2k rt k _e_ 
o z e 


^x^mx + ^y^my 



1-e 


-j2k z t 


(3.21) 
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-* 00 


As k 


-A! 77 


x y 


, k z -jk, and thus one can write, for large k. 


r d 

T, - 


-^o 


i ■ v*. 


< E e k o - k x > 7x ' Vy % 


(3.22) 


T u 
, I ^ 


^0 


^k ke 
o e 


<S k o - k y > 7y 


k x k y J mx 


(3.23) 


and 


T a 

I o 


~ J n 0 
ke 

o e 


k J + k J 
x mx y my 


(3.24) 


Comparing Equations (2.51), (2.52) and (2.53) with Equations 
(3.22), (3.23) and (3.24), it is readily seen that choosing 


(er+1) 

e e = 2“ (3.25) 


will make T a = "T, a=l,2,3, for large k. In other words, if one 
considers the grounded, homogeneous, dielectric half space D with e g 
specified by Equation (3.25), then the following integrals: 


mn 


z mn “ z mn 


(3.26) 


-J °» 
4^ " 


< T ! - 


.cos<t> m „ - J mw i s i rt4> m „ 
nx T mn my T mn 


d i 

- (t 2 -t 2 ) 


J n x' 5in +mn +J ny' cos *mn 


- j(k x* B n + VnJ „„ .. 

> e J dk x dk y 
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and 


m 


’m 


m 


(3.27) 


-jli 

4 ^“ 


CO 

// 

_oo 


T -T d 
1 3 1 3 




f + Vf> 


dk x dk y 


will converge rapidly since the functions (T - T d ), a=l,2,3, vanish 

identically for large k. One can compute r and P according to 

mn m 

Equations (3.26) and (3.27), and z d and v d via the alternative forms to 
n mn m 

be discussed in the next section, whereby z and v can be evaluated 

mn m 

efficiently as 




(3.28) 


and 


v 


m 




(3.29) 


C. ALTERNATIVE FORMS OF z d and v d 

mn m 

In this section, expressions of z d and v d are obtained directly 

r mn m 

from the vector potential that employs the homogeneous space Green's 
function instead of its plane wave expansion. Thereby the resulting 
expressions are easier to compute, especially when the two dipole modes, 
or the dipole mode and the impressed source are well separated and t is 
not too smal 1 . 
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Referring to Figure 3.2, the vector potential due to dipole mode 
m can be written in terms of the Green's function of the homogeneous 
half space D as follows: 


Sm 


-j k_R 

„ J e mns ^ e mm 

e e 


'mns 


'mm 


J m (x,y) dxdy, (3.30) 


for r = xx + yy and f'eS , where 

mn - 7 - , mn n 


R = R 
mns mns 


= |x(x m „ + x'cosa - y'sin* - x) + 
1 mn r mn J mn 


+ y(y mn + x'sin<|) mn + y'cos^H - y)| 


T mn 


mn 


and 


R . = |R . 
mm 1 mm 


= + z2t I 

mns 


Making use of Equation (3.7), it can be shown (by direct calculation) 
that the electric field at r+r' due to mode m is given by 


- . _ J^e 
E m (r+r') = f dxdy { 


e S, 


m 




’ R mns) R ms^ R mns^ 
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( 3 . 31 ) 


A A 


- (Jjx,y)Q,(R.) + (Jjx,y)*R mn jR.Q~(Rm .)) 
v m v mni ' ' nr ' mm' mm Z m" 


> 


9 ? ~j^Q^ 

*— n \ 


where Q^R) = (l+jkgR-k^R") e c / R 3 


( 3 . 32 ) 


2„2> " jk “ R 


Q 2 (R) = ( 3 +j 3 k g R-kgR ) e c / R 3 


( 3 . 33 ) 


R = |R|, R is the unit vector in the direction of R, 


kg — ( iW | 1 q £q Sg and 0 g — ^ Uq/ £g 


From the definitions of z d and v d , it is seen that 

mn m’ 


z mn = “ ; < *x , dy , 3n(x‘ ,y 1 ) • 


Jn e 

4 rt 7 { dxdy {( J m (x ’ y)Q l (R mns ) + 
e Sm 


A A 


+ (J m (x,y) *R mnc ) R rT ,„ f .Q 0 (R m „„ ) ) - (J m (x,y)Q, (R . ) 
m v ' inns' mns 2 mns' m ’•''u' mm' 


, + 


A A 


+ ^ d n/ x,y ^ * R mni^ R mni ^ R mni^ 


( 3 . 34 ) 


and 


v_ = 


/ dzJ. (x f ,y f ,z) 


jn e 

WT ' dXdy K^U.ylQjOVO + 

e S m 


A A 


+ • R fJ R fsV R fs>> - ' 3 m< x - y )0l< R fi> + 
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(3.35) 


where J.j (x^ ,y^ ,z) is defined in Equation (3.15), 

R fs = x(x f -x) + y(y f -y) + zz , R fs = |R^ S | 
and 

R fi = x(x f -x) + y(y f -y) + z(z+2t) , R fi = |R fi | 

It may be remarked that in computing z^ n and using Equations 
(3.34) and (3.35) one will encounter some numerical difficulty (namely, 
the loss of accuracy due to the differencing of two terms of almost 
equal values) when t is small. In such case, it is appropriate to 
expand the Qj and Q 2 functions as Taylor series, and to keep the first 
few leading terms. Expansions about z=-t may be most adequate. 
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IV. SUMMARY 


This report has presented expressions for: 

1. the self or mutual impedance between microstrip expansion and 
test modes 

2. the mutual impedance between a microstrip mode and a vertical 
current filament in the dielectric. 

These are the main quantities needed in an MM solution of the microstrip 
antenna. Methods for the efficient evaluation of the impedances are 
presented. Finally, expansion modes are described which would be useful 
in analyzing the microstrip antenna over a wide frequency band. 
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APPENDIX A 


DETERMINATION OF SPECTRAL FUNCTIONS 

Expanding Equations (2.2) and (2.3) in rectangular coordinates and 
making use of (2.4) and (2.5) yields the following: 

In region 1 (dielectric). 


3 l g2 

E xl (r) ■ - “ + 


3y v el jue i 3x3z Mml 

j j ~ , -j( k x x+k v y) 

4? IS Vel - ^7 VzlW e sir,k z ,(z + t)dk x dk 

-« J 1 


y ’ 

(A.l) 


_3 1 3 2 

E yl^ “ 3x ^el + jo)e 1 3y3x ^ml 


-j j ~ -j(k„x+k y) 

= 4^ II (Vel + ^ k v k zlW e sink zl (zn)dk x dk , 


y 

(A. 2) 
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E zl (?) = 


_ jO_ 
4tt 2 




_ J 

4ir2 


Hyi(F) - 


_ ii_ 

4tt 2 


H zl< ? ) ■ 


_ ~j 

4tt z 


3 Z 


[ 9^2 + £ r k n 1 

1 *• •• 


jwe, L 3z^ ~r 0 J Mill 
.2 . ,2 


// 


K k # 

— • Lift 

toe. 


-j(k x+k y) 

%,j e y cosk 2l (z+t)dk x dk y 


1 '3 2 


jojp Q 3x3z v el 3y v ml 


(A. 3) 


fj 


jk x k zl 




*el " k y'"ml 


~j(k x+k v) 

cosk , (z+t)e y dk dk , 

z x a y 


(A. 4) 


1 3 2 3 

j ajy Q 3y3z ^el " 3X ^ml 


-j(k x+k y) 

e y cosk zl (z+t)dk x dk y , 


(A. 5) 


d 


Jk y k zl 




\p . + k ^ 
r el x 1 


ml 




3 2 2 

322 * ) *■ 


el 


// 


ai v 


l/ 2 2 ~ 

(k x + V *el 


-j(k x+k y) 

e y sink zl (z+t)dk x dk y . 


(A. 6) 
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In region 2 (air) 


E x 2 ( ? ) - 

_ j 

4it 2 


E y2 (F) = 

- ~j 

4ir 2 


e z2 (?) - 
_ ri. 

4ir 2 


h „ 2 (F) = 

_ j 

4tt 2 


1_ 3 2 

3y ^e2 + jtoe 2 8x3z 'W 


II {k y J e2 


k x k z2 


toe. 


*m2> e 


■j(l< x x+k y y+k 22 z) 


dk x dk y , 


(A. 7) 


a_ l_ a 2 

3x ^e2 + jtoe 3y3z *W 


// tfc x * e2 


k y k z2 ~ , _-j(k x x+k y y+k z2 z) 


toe ^ m2 ^ e 

2 


dk x dk y . 


Jtoe 




it + k2 J 


II 


toe. 


%& e 


-J'(k x x+k y+k z2 z) 


dk x dk y , 


(A. 8) 


(A. 9) 


l a 2 


j(op o axaz ^e2 + 3y %n2 


f I 


k x k z2 


^o9 " k .A 


toji 0 y e2 y y m2 


-j(k x x + y*k z2 z) 


dk x dk y , 


(A. 10) 


41 


1 g 2 g 

Hy2(r) “ gyg 2 tg2 - 9x ^2 

i “ 

= 4 ^ II 

— 00 

(A. 11) 


k y k z2 


+ KK 


a)u 0 v e2 T "x y m2 


-j( k x x+k y y +k Z 2 z ) 


dk x dk y , 




-j( k x x+k y y +k z2 z ) 


v eZ 


dk x dk y 


(A. 12) 


Enforcing the boundary conditions at z=0 (dielectric-air interface) 
specified by Equations (2.8) and (2.9) leads to 


<\1 

X 

LlJ 

II 

»-H 

X 

LU 


at 

z = 0 

E yi ‘ E y2 


at 

z = 0 

H x 2 - H xl 

°sy 

at 

z = 0 

H y2 + H yl 

■ °sx 

at 

z = 0 


Use of Equations (A. 13), (A.l) and 


(A. 13) 

(A. 14) 

(A. 15) 
(A. 16) 

.7) gives 
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k x k Z 2 ~ /v 

+ “73 = [ k ,A 


y v e2 T oje 0 %2 ^ y ^ el <oe 0 e r 
Similarly, use of (A. 14), (A. 2) and (A. 8) gives 


k x k zlVl sin ‘ : zl t • (*.17) 


~ kyk z 2 ~ r ~ j ~ i 

k 4 > „ - ip o = fk ip , + k k . i|> . 1 sink -t 

x y e2 we. r m2 L x r el uie.e x zl y ml J zl 
u up 


(A. 18) 


Use of (A. 15), (A. 4), (A. 10) and (2.10) gives 


k x k z2 „ 


Jk x k zl „ 

(i)y Q *e2 • k y*m2 


ojy 0 ^el k y^ml 


cosk zl t 


~j J sy * 


(A. 19) 


Use of (A. 16), (A. 5), (A. 11) and (2.10) gives 


k y k z2 . 


Jk y k zl „ 

„u 0 *e2 + k x*m2 

+ 

m 0 *el + k x' k ml 


cosk zl t = -jJ cv . 


sx 


(A. 20) 


Adding (A. 17) multiplied by k y to (A. 18) multiplied by k x yields 


< k x + k y^e2 = + k y>+el sink zl t 


+e2 = sink zi t 'f' e i 


(A. 21) 
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Adding (A. 17) multiplied by k x to (A. 18) multiplied by (-ky) yields 


( k x + k P k z2 „ 


-jk 


we. 


Zl 2 2 ~ 

(K + K) t,i sink ,t 


W " we e r y ' 111 

or J 


-jk 


zl 


%Z ■ 7^T Slnk zl l *ml 


(A. 22) 


Adding (A. 19) multiplied by k x to (A. 20) multiplied by (-ky) yields 


k z2 k zl 

— <k* + k y )+ e2 -j ^ <k* * ky)cosk 2] t » el - j[k y O sx - k x J sy ] 


t k x ^ 


[k zl cosk zl t5; el + jk z2 ; e2 ] = k x j $y - k y j sx . 


(A. 23) 


Adding (A. 19) multiplied by ky to (A. 20) multiplied by k x yields 


(k x + V*,2 * l k x + k y>*ml cosk zl t = -i[ k y J sy + k x J sx> 


< k x + VK 2 + V cosk zi t ] = - J lVsy + Vsxl 


(A. 24) 


Using (A. 21) in (A. 23) leads to 


wy. 


"el 


, - - [k x J sy - k J ] 

(k x + k v )0, sy y 


2 2 
C x "y J "e 


(A. 25) 
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*e2 


coy q sink ^ t 

,71 2 " [y„ - k v J sx ] 

( k x * k y )D e x sy y s 


where 


D e = k z i cosk z i t + J‘ k Z 2 si^zi 11 


(A. 26) 


(A. 27) 


Using (A. 22) in (A. 24) leads to 


7m 1 


~ J e r lc z2 

+ \K 


f k x J sx + k y^sy] 


-k^sinlc^t 
W = ( k x + k >m 


t k x J sx + k y J sy^ 


where 


D m = e r k z2 cosk zl t + jk zl sink zl t 


(A. 28) 


(A. 29) 


(A. 30) 


45 



APPENDIX B 


MUTUAL IMPEDANCE BETWEEN TWO PATCH DIPOLE MODES 

In this appendix, the exact expressions for the mutual impedance 

(z ) between surface patch dipole modes 3^ and 3 n on a lossy grounded 

dielectric slab is presented. The mutual impedance (v ) between dipole 

mode J and the impressed source current 3^ is also presented. These 

calculations are essential in solving a microstrip antenna problem by 

the moment method, z „ and v form the elements of the impedance matrix 

mn m 

and voltage vector, respectively. 

Figure B.l shows two dipole modes, J m and J n , located on the 
surface of a grounded dielectric slab with parameters u Q and e r e Q . The 
ambient medium is free space with parameters and e Q . Mode m is 
centered with respect to the (x,y) coordinate system. The center of 
mode n coincides with the origin of the (x ' ,y ' ) system which is 
displaced from the center of mode m by a position vector ( x 0 »y Q )* The 
x'-axis is at an angle a with respect to the x-axis. 
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Figure B.l. Two expansion dipole modes on a grounded dielectric slab. 
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The mutual impedance between modes m and n is given by 


z = - ( E m (r +r')*J n (r')ds l 
mn ' m x o n x ' 


(B.l) 


where E m denotes the electric field excited by J m , and the integration 
is over the surface of mode n in the z=0 plane with 

AAA A, 

r' = x'x' + y'y' = x(x*cosa-y 'sina) + y(x'sinct+y 'cosa) (B.2) 


r = xx„ + yy„ 
o o JJ o 


(B.3) 


Without loss of generality, it is assumed that 


J m (r) = xJ(r) + yJ mw (r) 
m mx x J my x ' 


(B.4) 


A A 


where r = xx + yy and 


J n< ? '> = x ' J n + *V (F ') 


(B.5) 


It then follows from (2.17), (2.18), (2.24) and (2.26) that, at z=0. 


k x k z2 ~ , -j(k Y x +k v ) 


E (r +r') = 77 // [k * 0 + * Je x 0 y0 

mx v o 4 tt 2 jj l y y e2 uje. mi2 J 

- co J n 


• e-Jt-Cx' -y' dk dk 
x y 


(B.6) 
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-j - ~ Mz 2 ~ , ~j( k y X n +k v-0 

VV r '> = 4 ^ II t k x*e 2 ‘ «e„ Vl e ^ 


. e -j+(x',y',<») dk dk 

* y 


where 


k_„ = / k Q -k*-k y , Imk z2 < 0, Rek z2 > 0 


l z2 


4>(x * »y 1 »cx) = x‘ (k coso + k sina) +y*(-k sina + k cosa) , 
X y x y 


~ u>y 0 sink z it ^ ^ 

^e2 _ (k*+k*)D [ k x J my " k y^mx-l 


x y' e 


-k z isink z it 

%i2 " (k^+k^D [ k x J mx + k y J my^ 5 


x y' m 


.A 


2 2 2 
k zl * / c r k o- k x- k y 


D e = k zl cosk zl t + J k Z 2 sink zi t 


D m = e r k z2 C0Sk zl t + ^zl^^zl 1 • 


J my( k v» k J ■ II J mx (x ’^ )e 


j(k x x+k y y) 


tux' x* y 


dxdy 


(B.7) 

(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 
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and 


j(k x+k y) 

^my^x’^y) = ^ J my (x,y)e dXdy * (B.16) 

Also one may note that the tangential component of E m (r 0 +F') on the 
surface of mode n is given by 

E^(r +F* ) = xE (r +r‘) + yE (r +r‘) 
m v o ' mx v o' 3 my' o ' 

= x ' [E mx^o +f:, ^ COSa + E my( f o +f:, ) sincJ 

+ y,[_E mx (F o +F,)sina + E my (F o +F,)coSa] * (B.17) 


Combining (B.l), (B.5) and (B.17) yields 


h n w n 


mn 


■ - / / l[ E mx ( V ? ' ,cosat E my( F o +i: '> sin “) J nx l(x '- J '' ) 


-d n _w n 


+ t" E mx (F o +F,)sina + E my^ F o +F ' ) C 0 Sa ] d n y 1 ^ x * ,y ' ^ l d y' dx ' 


ny 


Next let 


(B.18) 


^nx' (k x’ k y’ a) = 5" "" J nx -( x ^ , ) e “ J * (x, ’ y, ’ a)d y ,d x' 

“hn “Wr 


'n - "n 


*)n ^n 


= I I V< x ‘>*'> 

-h n -w n 

-j[x‘ (k x cosa+k sina)+y' (-k x sina+k cosa)] 


• e 


dy'dx' 


(B.19) 
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and 


^ny'^x’^y*^ = /" J ny • Cx ' ,y 1 )e * J 4 ' ( x ' ‘ * a) dy ' dx ■ 

“^n ~ w n 

^n w n 

= / I J ny . (x* ,y * ) 


~^n “ w n 


-j[x' (k v cosa+k w sina)+y' (-k^sinct+k, cosa)] 
e x y X y dy'dx' _ (B-M) 


Use of (B.6), (B.7), (8.19) and (B.20) in (B.18) leads to 


Sin ~ 4ir 2 


II 


k x k z2 


(y e 2 + Oieo t2^ C0Sa 


* *** 
- (M 


k y k z2 


x y e2 a)G 0 Si2 


i?) sina 


^nx' ( k x ,k y ,a ) 


k x k z2 


k y k z2 


( k y%2 + u£g ^ sina + ( k x*e2 “ u>e 0 i 2 ^ C0Sa 


J ny' ( k x‘ k y ,a ^ 


"j( k x X o + Vo^., 
e J dk x dk y 


(B.21) 


z mn ™ Equation (B.21) represents the most general expression for 
the mutual impedance between two rectangular dipole modes on a planar 
grounded dielectric slab. 
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Referring to Figure B.2, the mutual impedance between dipole mode 
and impressed source is defined as 


v = / E_(r) *J. (r)dv 
m y m' ' l ' ' 


(B.22) 


Since source J-j can be reasonably modeled as a vertical filament of 
constant current inside the dielectric, it is represented by 


J.(r) = 2 l i 6 (x-x f ) 6 (y-y f ) , -~<x,y<» , -t<z<0 


(B.23) 


where I-j is the magnitude of the feed current and (xf,yf) denotes the 
filament location in the x-y plane, with respect to the center of 
mode m. 

Making use of (2.13), one can write 


_ 7 <V k y> ~ -^W 1 

mz^ r - 4ir 2 1 ’ ye £ ^ml e 

-oo or 


cosk . (z+t)dk dk 

c. X J 


(B.24) 


Inserting (B.23) and (B.24) into (B.22) and carrying out the volume 
integration gives 

-ih *, ( k x +k y) ~ -J'( k y x f +k v yf) 
v - =“// — , l. ’Ll e y sink_,tdk..dk.. . 


m 4tr 2 J 1 ojeek. ^ml 
-« o r zl 


'zl x y 


(B.25) 
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Figure B.2. A dipole mode 3 m and an impressed current source 3.. . 
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(B.26) 


In both (B.24) and (B.25) 


jer k z2 


♦ml = - (kJ+kjX tkxJ,nx + kyJ ” yl 
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APPENDIX C 


PATCH DIPOLE MODES 

In this appendix, the specific patch dipole modes and their 
corresponding Fourier transforms are presented. In using the moment 
method to analyze the performance of a microstrip patch antenna at 
higher order resonant frequencies, it is appropriate to employ the 
entire functions as the expansion dipole modes. On rectangular patches, 
the x polarized current modes will have the form 

* , rnir H* 

J mn (x,y) = x {sin 2L^ (*+L x ) cos 2L^ (y+L y )} > (C.l) 

where m=l,2,3 n=0,l,2,3 -L < x < L and -L < y < L with 

x * y y 

L x * 0 *. L^. It should note that the dipole modes J mn have double 

subscripts instead of single subscripts which appear in Chapters 2 and 

3. The double subscripts allow the modes to vary as simusoidal 

functions in both the x-and y-coordi nates. 

It is easy to show that the Fourier transforms of J assume the 

J mn 

following forms: 
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■vvv 




= /• 

-L 


' 3 mn (x '->'> e 

" L x 


j(k x x+k y y) 


dxdy 


“ x V k x> Y n< k y> 


where 


imr 


V k x> * < 


cosk x L x 


, m- 1,3,5,.. , 


mir 

*- 2L x k x 


jmir 

~ sl " k x L x 

L x 


, m=2,4,6,... 


Itlir 

^L x k x 


and 


j2cosk y L y 


, n-1,3,5,... 


1- 2L y k y 


V k y> = 1 


2s1nk y L y 


. I JHL_ 
k y I !- 2L y k y 


2 ~\ 


, n=0,2,4,, 


(C.2) 


(C.3) 


(C.4) 
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